NASA/TM— 2002-211703 



Real-Time In Situ Signal-To-Noise Ratio 
Estimation for the Assessment of 
Operational Communications Links 


Robert M. Manning 

Glenn Research Center, Cleveland, Ohio 


August 2002 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
data bases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at h ttp:lhviviu.s ti. nns a. gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/TM— 2002-211703 



Real-Time In Situ Signal-To-Noise Ratio 
Estimation for the Assessment of 
Operational Communications Links 


Robert M. Manning 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 


August 2002 



Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / /gltrs.grc. nasa.gov /CLTRS 



REAL-TIME in situ SIGNAL-TO-NOISE RATIO ESTIMATION 
FOR THE ASSESSMENT OF OPERATIONAL 
COMMUNICATIONS LINKS 

Robert M. Manning 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


INTRODUCTION 

The assessment of the propagation conditions that prevail on a communications 
link, whether it is a point-to-point terrestrial link or an earth-space link, is of vital 
importance for the optimal operation of the link. The significant indicator of the 
communications quality of the link is, of course, the signal-to-noise ratio (SNR) y . 

If the value of this quantity goes below a certain given threshold due, for example, to 
atmospheric conditions such as rain, the bit error rate on the link becomes unacceptable 
for reliable communications integrity. Knowledge of the dynamic behavior of y is thus 
essential for the optimal implementation of procedures to mitigate the further degradation 
of y . In many instances, a separate propagation receiver is employed to acquire an 
associated signal transmitted by a beacon on the communications satellite. Measurement 
of the fading conditions of this signal level is then extrapolated to that of the communica- 
tions link. In an effort to ease the operational and financial burden to terminal operation, 
methods (to be briefly discussed below) have been advanced to estimate y of the 
communications link by measuring various quantities of the operational communications 
link itself. 

Many methods have been advanced to estimate the SNR using an active 
modulated communications channels. For example, the output of the receiver matched 
filter can be sampled, i.e., the voltage level V s for an output symbol, and the value is 

compared to a pair of a priori determined voltage levels ±a, a > 0 a < . The 

statistical frequency of occurrence n F of values which fail to fall within the interval 
[-a, a] is calculated as well as the total number n T of samples considered. The ratio 
rip/rij is related to the symbol error probability by 



where E s is the symbol energy (i.e., energy per symbol), N 0 twice the noise power at the 
output of the matched filter and Q(- ■ •) is the well-known <2-function. This relationship 
is then solved for the SNR y = E s /N 0 by using the inverse function Q 1 (•••). Although 
this method has several shortcomings, e.g., limited dynamic range and sensitivity to 
automatic gain control variations and inter-symbol interference, it suffers from an 
irreconcilable defect; the inversion of the function ()(■■■), necessary to obtain the estimate 
7 of 7 , is mathematically correct only if dQ/dy * 0 . In the event that dQJdy — > 0, the 
problem of determining 7 becomes ill-posed 1 . That is, a small error in the estimate of 
n F /n T leads to a large error in the estimate 7 . In fact, the ill-posedness of this problem is 
the major source of the lack of dynamic range. 
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An entirely different approach can be imagined which avoids the ill-posedness 
of that above and thus tends to be more robust in the presence of measurement errors 
inherent in the sampling of the matched filter output. In the BPSK case, the sampled 
voltage V s is related to the bi-phase signal amplitude ±A and the corresponding in-phase 
noise component A c by 

V S =*(±A + AI C ) 

where K is a constant proportionality coefficient incorporating receiver gain factors, etc. 
Remembering that the goal here is to obtain an expression for the SNR y = E s /N 0 

= A 2 /2o 2 n where a 2 N = (n 2 c } for zero mean, white Gaussian noise, the method endeavors 
to obtain this ratio solely from the measured values of V s . Thus, to separate the noise 
term, one would try to form the average (V s ) using the fact that (/V„) = 0 . However, the 
random bipolar nature of the signal amplitude A also yields a zero average giving 
(V s ) = 0 . The technique that is then adopted in this approach is to form the absolute 
value |P S .| of each sample and then forming the ensemble average giving 

{\V S \) = (\K(±A + N,f) 

~ KA 

so long as the condition A » o N prevails. Additionally, the sampled values V s are used 
to compute the variance 

{V s 2 ) = K 2 (A 2 +< Ti) 

Hence, using the former expression to rid to the A 2 term to give 

^ = <^ 2 >-(N} 2 

therefore allowing one to write 

.. ^ (M. 

2a » 2(N ! }-(M) 2 ) 

solely in terms of the voltage samples V s . 

The major drawback of this method is the formation of the absolute value of 
the random quantity V s ; such an operation can drastically change the statistical 
characteristics of the random variable. Even though this was done to rid of the bipolar 
nature of the communications signal amplitude, it is mathematically faulty and, 
moreover, needless. One can, and in fact should, incorporate the bipolar nature of the 
signal amplitude into a rigorous statistical analysis; this characteristic of the signal is 
just as important as its other aspects. In addition, this approach requires one to use a 
predetermined bit-stream format within the communications data composed of a series 
of l’s thus necessitating a synchronization with, for example, a preamble within the 
modulation format. This results in a further complication of its implementation and, as 
mentioned above, is unnecessary. 

It is the purpose of this work to formulate the rigorous statistical basis for the 
correct estimation of BPSK signal SNR from what is known about its behavior at 
the input and output of the receiver demodulator. Instead of employing tacit and 
unwarranted assumptions concerning the nature of the communication signal for 
analytical simplification, a complete consistent statistical description of a BPSK signal 
will be provided to which the well-known techniques of maximum likelihood estimation 
theory can be applied. By employing, rather than neglecting, all the subtitles of the 
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statistics describing the BPSK signal, an unbiased estimation procedure will be derived 
that makes simple use of its inherent phase characteristics at the demodulator. In what is 
to follow, a preliminary review of BPSK signal representation will be given which will 
lay the foundation for the statistical connection between Gaussian noise and SNR. Once 
an appropriate probabilistic description is obtained that establishes a rigorous contact 
between SNR and the measured phase error of the BPSK signal entering the receiver 
demodulator, the methods of maximum likelihood estimation theory will be used to 
obtain analytical expressions for biased and unbiased estimates of SNR from easily 
measured phase errors. Finally, the straightforward modifications needed at the 
demodulator to implement the required phase measurements will be given. It should 
be noted that the resulting SNR estimation technique is also applicable for a QPSK 
demodulator simply by applying it to one of the BPSK arms with appropriate 
modifications for the SNR expression. 

PRELIMINARIES OF SIGNAL REPRESENTATION 

Consider the simple BPSK demodulator shown in Figure 1 . The phase-modulated 
signal s t (t) at the input is defined by 


5 i .(t) = A(t)cos(m o t + 0, : (t)), , i = 1-2 ; M = 2 (1) 

where co 0 is the angular frequency of the carrier wave, A{i) is its time varying amplitude 
and 6 t (t) is its bi-phase state. In terms of quadratures, the noise-free BPSK signal is 
given by 

5 ; (t) = A(f)[cos(m,/)cos(6I(f)) -sin(w 0 t)sin(6>(f))] 

= A(t)cos(ft) 0 r)cos(0,(t)) 

= A c (t)cos(m 0 t), A c (t) = A(t)cos(0 f (t)) (2) 

where 

A c (0 s A(t)cos(0,(t)) = ±A(t) (3) 

showing the obvious fact that the noise free signal only has an in-phase component A r {t) 
since the bi-phase states are #,(/) = 0°, 0 2 (t) = 180° thus relegating the sine factors to zero. 


2cos(a>ot) 


e(t)= Sj ( f) + n(t) 




In-Phase Output 

► 

A c (t) + N c (t) 


s { ( t ) = A(t) cos (co 0 t+ Oi (0) A c (0 = A(f) cos (Of (t)) = ±A(t ) 

n( t ) = N(t) cos( ffip t + 0 N ( t)) N c (f) = N(t) cos(0 N ( t )) 


FIGURE 1. 

BPSK demodulator showing relationship between signal parameters 
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Similarly, the noise n(t) associated with the input is given by 

n(t) - N(t ) cos [o) 0 t + 0 N (t)) (4) 

where N(t) is the time vaiying amplitude of the noise and d N (t) is the associated phase. 
In terms of quadratures, the noise is given by 

n(t) - A^(t)[cos((u o r)cos(0 v (t)) -sin((u o r)sin(0 A ,(r))] 

= N c (t)cos(co 0 t ) - N s (t)sin(o) 0 t) (5) 

where 

N c (t)^N(t)cos(d N (t)) (6) 

and 

N s (t) = N(t)sm(d N (t)) (7) 

showing that, unlike the noise-free BPSK signal component, the noise is characterized by 

both in-phase and quadrature terms. The composite signal e(t ) at the input of the 
demodulator is given by 

e(t) = s^t) + n(t) 

= [A c (t) + NXt)]cos(co 0 t) + N s (t)sin(co 0 t) (8) 

where use has been made of Eqs.(2) and (5). Writing the composite signal e(t ) in terms 
of its in-phase and quadrature components and employing Eq.(8) yields 



e(l) = £(/|cos(«/) - C. sin(® 0 () 


(9) 

where 

£4<)- 4(0 +«,(<) 


GO) 

and 

£.(0 = 4(0 


(ID 

Writing Eq.(9) in the standard form, in which the signal and noise components are written 

in Eqs.(l) and (4), finally gives 




e(t) = Ecos(d) 0 t + 0 E {t)) 


(12) 

where 





E(t ) = ^E;(t) + E;(t), e E (t) = tan 1 

P 

1 

(13) 

Thus, 

E c {t ) = E{t)cos(6 E (t) ) 


(14) 

and 

CD 

ii 

bf 


(15) 


The geometrical depiction of the various quadrature components dealt with here in the 
case of the bi-phase signal ^(f) = A c (t ) = A(f)cos(0j(t)) = +A(r) is shown in Figure 2. In 
this figure, it is easily seen how the addition of the noise n(t) to the signal ,v, (t) yields a 
total composite signal ^(t)with a phase of d E (t ) with respect to the in-phase axis. The 
phase 0 E (t ) becomes, in this bi-phase signaling case, the phase error of the received 
signal .v,(r), which represents the input to the demodulator. 

Given the preceding development of the relationships of the various signal 
components that enter into the demodulation of a BPSK signal, the problem to be 
addressed here is the determination of the associated signal-to-noise ratio 
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E c 


FIGURE 2. 

Relationship of the various signal and noise quadrature components of the input of the 

BPSK demodulator 



As shown in the next two sections, this problem can be placed on a 
mathematically tractable basis and its correct solution rests on the statistical connection 
between the easily measurable parameters, not of the demodulator output, but of the 
demodulator input, i.e., E(t) and/or d N (t), and the signal-to-noise ratio y. The fact that 
the demodulator input is significant is that it contains both the in-phase and quadrature 
components of the noise. Additionally, the bi-phase nature of the signal modulation, i.e., 
±A(f), is a major impediment in the application of the methods mentioned above. This 
aspect of the signal must be properly accounted for in a legitimate treatment of the 
problem. Such a program begins with the use of what is known about the statistics of 
the noise process. 
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NOISE STATISTICS AND THE SIGNAL-TO-NOISE RATIO 

The noise components N c (t) and N s (t) are random quantities whose statistics are 
usually taken to be governed by a Gaussian random process. Additionally, it is known 
that such quadrature components are statistically independent. Letting a 2 N (t) be the 
variance of each of these noise components, one has for their joint probability density 


p (nmnM) = 


i 


exp 


Nm+im 

2 oUt) 


(16) 


.2^(0 

Using Eqs.(10),(ll), (14) and (15), one has in terms of the composite signal parameters 
E(t) and 0 E (t), 

N c (t) = E(t)cos{d E (t))-A c (t) (17) 

and 

N,(t) = E{t)sm{0 E (t)) (18) 

Transforming the probability density function of N c and /V s (dropping the time 
dependence for notational clarity) given by Eq.(16) into one that is a function of E and 
0 E {t) involves the transformation 

b(y-'y,)l 


p(E,e s ) = p(N',N,) 


which requires the Jacobian 




d{E,6 E ) 


dN,. 


d{E,e E ) 
dN, 


(19) 


dE 

dN, 


dO E 

dN, 


dE d0. 


= E 


( 20 ) 


where Eqs.(17) and (18) where used to obtain the indicated result. Thus, using Eqs.( 1 6)- 
(20), one has in terms of the composite signal parameters 


p{e,o e ) = 


\ 


2na 


exp 


N J 

E 


(EcosO e - A c )" + E 2 sin" 0 E 


2a: 


2 k a 


exp 


N J 


|2 X 


V 


2a 


exp 


n J 


E " -2A c EcosO e 
2 ai 


( 21 ) 


The signal-to-noise ratio y characterizing the composite signal input to the demodulator 
is given for the BPSK case as y(t) = A(t)/a 2 N (t). Writing Eq. (21) in terms of y and 
remembering Eq.(3) gives 


p(E,0 E \y,0i,a N ) = 


2na 2 N j 


exp 


7 


cos‘ 6, 


•exp 


E 2 -2Ea N - s fy cosO t cosO E 

2^ 


( 22 ) 


where the density function p(E,0 i: ) is now written as a conditional probability density 
/?(.E,0 £ |y,0 1 ,<7 A ,) governing the values of E and 0 E conditioned on the values of y, 0,. 
and a N . This is done with the hope of being able to obtain an expression, using Eq.(22), 
connecting the measurable values of E and/or 0 E at the output of the demodulator only 
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to the associated value for y thus allowing one to statistically estimate y from such 
easily measured values. 

Having secured the above relationship governing the probability density of the 
envelope and phase of the composite signal, i.e., E and d E , and the prevailing values of 
the signal-to-noise ratio y, the variance of the noise power <j 2 n , and the phase state of the 
transmitted signal 9 t , one can obtain two relationships involving measurable statistics of 
E, or d E , and 7 . For example, one can obtain the statistics, conditioned on the value 
for 7 , that describe the composite signal envelope over all possible signal phase states 
simply by summing Eq.(22) over the two possible values for e i (t) and integrating over 
all possible values of the phase d E , viz 2 .. 


J p[E,9 E \y ,9 t , °n = 

1=1 


( E A 
\2no n j 


ex P| -yJ ex P| 


f E 2 ^ 


V ~°N J 


£ 


exp 

+^ C ose E 

+ exp 

A £ cos9 1] 

l 

L 


L °N JJ 


dO, 


E f 7 

— ex Pl — T l ex P 


V ^ N 


( 


2 < t ;, jj 0 


f 2 E^ 


_2 

\°N J 


Y 


ex P| ~ | exp 


f 


V 

E 2 


1 


K f 

2 cosh 


4ye 


V °N 


cos 6 f 


de F 


2o: 


AM. 1 

i ^ J 


= p(E\y,o N ) 


(23) 


which is a form of the well-known generalized Rayleigh distribution for the bi-phase 
signal envelope. Although this expression is useful for many purposes, it is not sufficient 
for the estimation of the quantity 7 from measurable statistics involving E(t) since, as 
can be seen from Eq.(23), it explicitly involves the parameter a 2 N . Thus, using this 
approach, one must know, a priori , the prevailing values of a 2 N . Further consideration 
shows that this must be the case since the noise n(t) has two degrees of freedom which 
contributes to the composite signal envelope, i.e., N(t) and d N (t), as seen from Figure 2. 
Integrating over all possible values of the phase characterizing the demodulator output 
leaves the phase of the noise to be determined by knowledge of o N . Thus, if one desires 
to obtain an estimate of 7 from measurements of E or some statistic related to E, one 
must also need an a priori estimate of o N thus rendering useless a straightforward 
estimate of 7 from easily measurable signal parameters. 
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However, if one considers statistics that describe the composite signal phase error 
9 E (t) by summing Eq.(22) over the two possible values for (9,(7) and integrating over all 
possible values of the signal envelope E(t), one obtains 3 

J 0 S />( £ '0ir|r.0,.<7»)<! , E = 



= p(« E |r) (24) 

where the third integral results from the change of variables r = E/o N and erf (• • •) in the 
fourth line is the ‘error function’ or probability integral defined by 

erf(jc) = -^=Jexp {-t 2 )dt 

This formulation does not involve any a priori information other than the signal-to-noise 
ratio y . Thus, Eq.(24) gives a relationship involving the conditional probability density 
of the phase of the demodulator output given a value for the signal-to-noise ratio y of the 
composite signal, with no other a priori information necessary. This expression is indeed 
a candidate for the basis of estimating y from knowledge of values for 6, . The proper 
mathematical foundation for this procedure will now be given. 


THE MAXIMUM LIKELIHOOD ESTIMATION OF SIGNAL-TO-NOISE RATIO 
FROM PHASE MEASUREMENTS OF THE BPSK DEMODULATOR INPUT 


Equation (24), giving the conditional probability density governing the values of 
6 e given a prevailing value for y, can formally be written in the opposite sense, i.e., a 
conditional probability density governing the values of y given a prevailing value for 0 E , 
viz., 



Although not explicitly shown in Eq.(25) for notational simplicity, all signal components 


that enter into this expression are still functions of time, 0 E = 0 E {t), etc. Assuming that a 
time series of statistically independent phase measurements 0 E [t J , d E (t 2 ), •••, d E (t k ) 
can be obtained from the composite signal input to the demodulator that subtend a time 
interval At = t k - f, whose value is small enough such that y(t) and a(t) can be taken to 
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be constant in At, one can use Eq.(25) to obtain a probability density function 
conditioned on the series of phase measurements given by 

k 

/>(r|<MO. ■ ■ -AM) = n^Mo)) (20 

;=i 

From this, one can form the likelihood functional defined by 

L (r|M0’ M* 2 )> ••• A(O) = ln {p(y|MO’ MM -"MO)} 

k 

=5>{/>(r|<M'4 (27) 

7=1 

According to the method of maximum likelihood, the corresponding estimate y * of y is 
found from this functional as that value of y for which L(- ■ ■) is a maximum. Hence, the 
estimate y * of y which prevails over the set of phase measurements Mf ), 0 E (t 2 ), 
0 E (t k ), of a BPSK modulated signal is given by 

MV|MO’ MM '"’MO) _ A 

- U (zo) 

dy 

7 = 7 * 

It is important to remember that the bipolar nature of the modulated signal as well as the 
usual assumptions of Gaussian noise are already convolved in the probabilistic 
description of Eq.(25). Also, unlike the general case where y and a 2 N are functions of 
time, they can now be taken as constant and independent of time during sufficiently small 
measurement intervals At. A quantitative measure for ‘sufficiently small’ will be given 
below. 

It is now necessary to find the root y * of Eq.(28). Using Eqs.(25)-(27) in 
Eq.(28), and, for analytical tractability, neglecting the last term of Eq.(25) (an 
approximation which holds for large values of y ), and performing the required 
differentiation yields 

sin 2 (MO) + 


Completing the differentiation of the last term results in a rather unwieldy expression. 
Using the fact that the erf(-- ) function tends to a constant for large values of the 
argument, i.e., for large y, one can neglect this term (consistent with the neglect of the 
second term of Eq.(25)) and find that in this case, the optimal estimate for y of a bi- 
phase signal with additive Gaussian noise based on a series of k phase measurements 
MO °f the composite signal input is given by 

i * V 

-Ysin 2 (M0) > y»\ (30) 

;=i 

Because of the approximations involved with neglecting the second term of Eq.(25) and 
the third term of Eq.(29), this maximum likelihood estimate is biased toward large values 
of y . This bias must now be removed for small values of 7 . To this end, one uses the 
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original expression of Eq.(25) retaining the second term since this will be appreciable for 
small y . However, to maintain analytical flexibility, the erf(---) function of Eq.(25) will 
continued to be neglected. Finally, the analysis becomes more amenable if the inverse 
function \/y * given by Eq.(30) is used. Given these considerations, one has 4 for the 
unbiased estimate y of y * over all possible values of the phase 9 E (t) 



Thus, the unbiased estimate y of y is related to the biased estimate y *, given by 
Eq.(30), through the non-linear relationship of Eq.(31). (It should be mentioned that the 
ensemble average of sin 2 (0 £ ) is calculated over the phase interval -n/2 <9 E < n/2 since 
p{y\0 E ) is periodic in d E with period n. (Hence, phase values that differ by ±n 
essentially correspond to the same SNR values and the averaging interval chosen in 
Eq.(31) implements the separate identification of these intervals.) Since there are no 
other parameters of the problem that enter into Eq.(31), this relationship is universal for 
BPSK modulation and is plotted in Figure 3. From Eq.(31 ) one has the following 
limiting behavior: 



Thus, for large y , the two estimates converge for the reasons explained above. In the 
case of small y, the biased estimate approaches 2 since in this case, the phase errors 
9 E {t t ) are randomly distributed within their range -n/2 < 9 i: (t i ) < n/2. Note that Eq.(31) 
does not monotonically approach the value of 0.5 as y 0 ; it overshoots the 0.55 level 
before it turns toward the 0.5 limit. This defect is due to the approximation made above 
in the neglect of the erf(- • •) function in Eq.(29). A more careful analysis may correct this 
shortcoming. 

The implementation of this process is straightforward. Considering now the 
composite signal input to the demodulator in terms of Eqs.(9), and (14)-(15), instead of 
the previous representation of Eq.(8), and keeping in mind the goal of securing the values 
of the phase error 0 E {t), one arrives at the required modifications shown in Figure 4. As 
k sampled values of 9,{t) are obtained (just how to determine a value for k will be 
discussed below), one then implements Eq.(30) to obtain the biased estimate y * of the 
signal-to-noise ratio y . Once this has been obtained, Eq.(31) is employed and solved for 
the corresponding unbiased estimate y ; this latter procedure can be effected by use of a 
‘look-up’ table that represents the universal graphical behavior of Figure 3. 
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7 

FIGURE 3. 


Characteristic curve relating the biased estimate y *, calculated from phase error 
measurements d E {t t ) according to Eq.(31), to the unbiased estimate y . 


DETERMINATION OF PHASE SAMPLING LENGTH TO ASSURE 
CONSTANCY OF y AND a 2 N DURING THE INTERVAL 

It is imperative that the prevailing values for y and a 2 N do not evolve during the 
time period in which the k samples of composite signal phase are formed to obtain 
the biased estimate y*. The fundamental constraint which must be addressed during the 
measurement process is the fact that y does not have spectral components that exceed 20 
Hz for an atmospheric channel. (This is the upper limit for atmospheric scintillation at 
communications satellite frequencies.) Thus, if the data rate of the communications link 
is R d bps, and if a measurement of the phase 0 £ (f ( .) occurs for k consecutive bit 
intervals, one must satisfy the Nyquist sampling constraint 

— > 40 Hz 
k 

This inequality serves to bound the number of phase samples. 


NAS A/TM— 2002-2 1 1 703 


11 




2 cos( coqO 



0£(O 


= N s (t) 

FIGURE 4 

Modification to BPSK demodulator to obtain phase error in composite signal. 


EXTENSION TO THE GENERAL M-ARY PHASE MODULATION CASE 

The simplest application of the foregoing to QPSK and other M-ary phase 
modulated cases involves using the method shown in Figure 4 in conjunction with 
Eq.(30) with just one of the arms of the demodulator. In the interest of maintaining 
mathematical rigor, one could return to the analysis given above and employ, e.g„ M=4 
in Eq.(l) for the QPSK case. Thus, Eq.(24) would then employ a sum over i = 4 possible 
phase states. This, however, would complicate the SNR estimate analogous to that of 
Eq.(30) for the BPSK case. 
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